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ABSTRACT 


This paper presents an accurate Lithium-ion battery model representation in 
Matlab/Simulink. The Tremblay's battery model was used as a BES model 
platform, where the determination of the model parameters was obtained 
based on heuristic optimization approach. This approach is simple but more 
accurate compared to the conventional method. In the classical method, it 


requires the user to manually select the battery model parameters from 

relevant points on the manufacturer discharge curves. However, this way of 
Keywords: battery parameters extraction normally exposed to the human error and 
would easily result in an inaccurate selection of battery parameters for the 
BES simulation studies. Therefore, an easy and accurate approach using 
heuristic optimization for determining battery model parameters was 
introduced. The simulation studies utilized three different optimization 
algorithms for comparison purposes, i.e. 1) Particle Swarm Optimization 
(PSO), 2) Gravitational Search Algorithm (GSA), and 3) Genetic Algorithm 
(GA). The performance of BES model discharge accuracy with respect to the 
test data from three different algorithms was compared and the results 
showed that the GA approach gives the best results in terms of accuracy and 
execution time. Finally, the validated results of GA-optimized battery model 
showed the accuracy of 98% compared to the conventional approach. 
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1. INTRODUCTION 

Battery energy storage (BES) system is an important element in the renewable energy (RE) system 
and electric vehicles (EV) applications [1, 2]. In solar photovoltaic (PV) system application, the integration 
of BES with the PV system is one proven method in mitigating the output power fluctuations of PV sources 
[3]. Currently, there are various types of batteries commercially available for BES such as Lead-Acid (LA), 
Lithium-ion (Li-ion), Nickel Metal Hydride (NiMH), Nickel Cadmium (NiCd), Sodium Sulphur (NaS) and 
many more [4, 5]. Li-ion batteries have outstanding applications especially in the operation of plug-in hybrid 
electric vehicles (PHEVs) and power grid applications. This is due to its characteristics of high energy 
density, high charge and discharge ability, low self-discharge loss, longest cycle life, highest efficiency, and 
no memory effects compared to the all other battery technologies [4, 6]. 

Recently, a simulation study is a frequently used technique to design and develop a high efficiency 
and robust control strategy of the BES system, particularly in electric vehicles and power system 
applications. The simulation technique can reduce the cost of commercialization of new technology as it can 
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avoid the unnecessary procedures in testing so as avoiding the purchase of expensive measuring instruments. 
Simulation studies of BES control system have been extensively carried out in the past [7-13]. In this regard, 
several battery models were proposed to further evaluate and develop the BES control system. An accurate 
battery model is needed as it has direct influence to the state-of-charge which can affect the robustness of the 
BES control system [12-14]. 

In [15], a brief overview of several battery models for various applications are presented. The 
battery models are classified into several categories: simple model, Thevenin-based model, impedance-based 
model, runtime-based model, combined electrical circuit-based model, and generic-based models. A generic 
battery model is accurate and appropriate for all types of battery energy storages. In [16], an improved and 
easy to use generic battery model is presented for electric vehicle applications. However, the parameters of 
the battery model need to be manually obtained from battery manufacturer discharge curve. This results in 
poor model performance when implemented in the simulation software. The accuracy of the battery 
parameters extraction has direct effects on the performance of the battery model. To accurately estimate the 
model parameters, one effective way is by performing parameter extraction from manufacturer discharge 
curves using optimization approach. A Quantum-behaved Particle Swarm Optimization (QPSO) and Particle 
Swarm Optimization (PSO) are used to obtaining the Li-ion battery parameters from the manufacturer 
discharge curve for the electric vehicle application have been presented in [17] and [18], respectively. 
Whereas, in [19], a parallel Jaya algorithm is applied to estimate the Li-ion battery parameters and the 
simulation results showed good performance of the developed algorithm. The specific optimization 
algorithms used in [17-19] however gives challenges to the user with limited access to the optimization 
codes. An easy to use guideline for such parametric optimization strategy should be developed considering 
widely available optimization approach. This paper presents parameters estimation strategy of the BES model 
using widely used heuristic optimization approaches such as Particle Swarm Optimization (PSO) and Genetic 
Algorithm (GA). In addition, Gravitational Search Algorithm (GSA) is also used for comparison purpose and 
diversify the results of optimization. The generic battery model [16] of Li-ion type is considered in which its 
parameters are extracted from the manufacturer discharge curves. Finally, based on the proposed heuristic 
optimization approaches, the comparison of Li-ion battery model discharge performance is presented. 


2. RESEARCH METHOD 

The present work is carried out from the development of accurate Li-ion battery model. Then, the 
parameters of the developed Li-ion model are optimized using heuristic optimization methods considering 
the PSO, GSA, and GA respectively. Then, the optimal Li-ion battery model is validated by comparing the 
characteristic curves of the optimal Li-ion battery model developed, with the test data from the manufacturer. 


2.1. Development of li-ion battery model 

The Li-ion battery model is implemented according to the Tremblay’s generic battery model 
presented in [16]. This model is high efficiency and can give good performance in studies related to electric 
vehicles and renewable energy technologies. The model can be represented by using equations to describe the 
electrochemical behavior of the battery in terms of state-of-charge (SOC), terminal voltage and internal 
resistance. The formulated equations of the dynamic model are described as follows: 








Veat = Ebat =i Rint! Bat (1) 
SOC = 100 (Hes) (2) 
Epataise = Eo — |K (-%-i*)| - |k (Sit) | + Ae (3) 
Epatcharg = Eo — |K (= i*)| - [k (Sit) | + ae (4) 


where Vga is the battery voltage, Rin, is the battery internal resistance, Igar is the battery current, Q is the cell 
capacity, Egar is the battery electromotive force, Exardisc and Egarcharg are the battery electromotive force 
during charge and discharge, Eo is the battery open-circuit voltage, K is the polarisation constant/polarisation 
resistance. The term it can be obtained by integration the battery current, i.e. f Isar dt, which is the actual 
battery current. Whereas, i“ is the filtered current. The exponential zone of battery discharge curves are 
represented by A, which is the exponential zone voltage, and B represents the exponential zone time constant 
inverse. For the battery model, there are several specific assumptions and limitations such as there is no self- 
discharge, the nominal capacity and internal resistance are constant, and there are no environmental 
considerations. From the developed battery model, BES system model can be constructed by series/parallel 
combination of single battery cell parameters. Table 1 illustrated the transformation of battery cell 
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parameters, where the number of cells in series (ns) and parallel (np) determines the total output terminal 
voltage and capacity or total size of a BES system, respectively. 


Table 1. Battery parameters transformation 


Parameters (unit) Value 
End of nominal zone voltage (V) Vnom X Ns 
End of nominal zone capacity (Ah) Qnom X Np 
Rated capacity (Ah) Qrated X Np 
Internal resistance (Q) Tint X Np/ Ns 
Maximum capacity (Ah) Qmax X Np 
Fully charged voltage (V) Viull X ns 
Nominal discharge current (A) LTnomdise X Np 
End of exponential zone voltage (V) Vexp X Ns 
End of exponential zone capacity (Ah) Qexp X Np 


2.2. Battery parameters extraction method 

The battery parameters can be approximated by using the manufacturer’s data by following the 
procedures given in [16]. Figure 1 illustrates the battery parameters extraction procedure from the typical 
discharge curves from the manufacturer. As illustrated in the figure, three important points namely fully 
charged voltage (Vj), end of exponential zone (Ve, Qe) and end of nominal zone (Vnom, Qnom) are 
manually obtained from the figure. From these points, the parameters of Eo, A, B, and K can be determined. 
The parameters accuracy from this approach is depends on the accuracy of the points marked on the curves. 
In practice, it is difficult to identify a correct point on the curves just by using visual analysis. This is because 
the raw data of the manufacturer discharge curves normally not supplied to the user. Thus, an intelligent 
approach is a promising way to overcome this issue. By using an optimization approach, a good model 
design can be achieved according to application needs and can avoid human error. Furthermore, more time 
can be saved during the parameters extraction procedure. 
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Figure 1. Battery parameters extraction procedure based on the discharge curve of 
panasonic Li-ion CGR 18650AF battery [20] 


3. OPTIMIZATION OF BATTERY PARAMETERS 

In the optimization process, six optimal parameters needed to be obtained which are Van, Vexp, Qexp, 
Vnom, Qnom and Rin. The objective of the optimization problem is to minimize the deviation of the developed 
battery model discharge curves with the real battery discharge curve data obtained from the manufacturer 
datasheet. In this regards, the objective function is determined according to (5), where the vector x and k 
represent the battery parameters and the index of the sampled data, respectively. While, Voaimoaa) and 
Vhat manu) represent the voltage data of the developed battery model and the voltage data of real Li-ion battery 
from the manufacturer. The optimization problem is solved by three different algorithms which are PSO, 
GSA, and GA for comparison purposes. The justification for selecting these three algorithms is due to their 
wide applicability and readily to be implemented in Matlab. These algorithms also have been proven to 
deliver good optimization accuracy, particularly in battery storage application in renewable energy 
sources [7, 12, 13, 21-23]. 
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OF (x) = min | f (Vrat moa (K) — Voatmanu(k)) at! (5) 


3.1. Particle swarm optimization (PSO) 

The PSO algorithm has been invented by Kennedy and Elbert in 1995 [24]. The general process of 
the PSO algorithm is as illustrated in Figure 2(a). PSO algorithm solves the problem in (5) by having the 
random population of particle solutions (battery parameters). PSO algorithm optimizes the problem by 
iteratively moving these particles around in the search-space according to simple mathematical formula as 
described in (6), (7), (8), where X f and X di are the previous and updated position of particle j; V;* and y 
are the previous and updated velocity of particle j; Pbestj and Goes: are the best position of particle j and the 
best position of the entire particle solutions; c; and c2 are cognitive and social learning rates; r; and rz are the 
random number between 0 to 1; ©, @min, and @max are represented the total weight, minimum weight and 
maximum weight; it; and itmax are number of current iteration and maximum iteration, respectively. For a 
robust and fast convergence PSO algorithm process, parameters of C7, C2, @min, ANd Wmax are set to 2, 2, 0.4 
and 0.9 as described in [12]. 
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Figure 2. Optimization procedures using PSO, GSA and GA 


3.2. Gravitational search algorithm (GSA) 

GSA is an optimization algorithm based on Newton’s famous law of gravity and mass interactions 
[25]. The simplified process of GSA algorithm is illustrated in Figure 2(b), where the process is started with 
the initializing the random agent population. In GSA, the updated position and velocity of the agent is 
calculated using (9) and (10), where 7; is the random number between 0 to 1 and a; is an acceleration of the 


agent at the current iteration, respectively. The acceleration of the agent is obtained using (11), where Ff and 


M hm are the total force that acts on agent and the inertial mass of the agent, 
k+1 — k+1 
a (9) 
+h, — 
xp =rV +4; (10) 
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a; = Ff / Mint (11) 
3.3. Genetic algorithm (GA) 

GA is an algorithm inspired by the process of natural evolution. GA algorithm has been successfully 
applied to a wide range of real-world problem of significant complexity. Starting with a randomly generated 
population (chromosomes), GA algorithm uses three main types of rules (selection, crossover, and mutation) 
at each step to produce a successor population for the next generation as illustrated in Figure 2(c). During 
selection rules, the parent chromosome that contributes to the population is selected for the next generation 
process. The selected parent chromosomes are re-combined to produce child chromosomes. These processes 


are iterated until the satisfactory fitness level is reached. 


3.4. Validation of an optimal battery model 

The optimal Li-ion battery model in the simulation is validated based on the typical discharge 
characteristics curve of Panasonic Li-ion CGR18650AF from the manufacturer following the procedures in 
[16]. During validation, the constant battery current are set at 0.2C (0.43 A), 1C (2.15 A) and 2C (4.3 A). The 
results are as shown in the results and discussion section. 


3.5. Simulation set-up for battery model parameter optimization 

Figure 3 illustrates the overall simulation set-up diagram for optimization of battery parameters. 
During the optimization of the battery parameters process, the developed battery model in Matlab/Simulink is 
linked to the PSO/GSA/GA algorithm in the Matlab/M-file as presented in Figure 3. Meanwhile, the typical 
discharge data from the manufacturer data is saved in Matlab/MAT-file. Before the optimization process is 
started, the population number, dimension number, and iterations are initially set to 25, 6 and 100 for all 
algorithms, respectively. The optimization process by PSO, GSA and GA algorithm is illustrated in Figure 
2(a)-(c). The processes are started with the randomly set of the initial population for each battery control 
parameters. The generated random parameters are evaluated by using a fitness function in (5). The process is 
repeatedly executed until the optimal parameters of battery are determined. 
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Figure 3. Simulation set-up of the battery parameters optimization 


4. RESULTS AND DISCUSSION 

The development process of the Li-ion battery model is started in the Matlab/Simulink. Then, the 
process is continued with the battery parameters optimization. Finally, the validation of the optimal battery 
model using typical discharge curves of a real battery from the manufacturer is carried out. The complete 
results of the battery’s parameters optimization are presented in the following sub-section. 


4.1. Optimization of battery model parameters 

Table 2 illustrates the comparison of the obtained battery parameters for the cases of manual and 
optimization (PSO, GSA, and GA) approaches. Besides that, the performance and elapsed time of 
optimization approaches also presented. From Table 2, the GA approach converged faster compared to the 
other approaches with the elapsed time of 1680 s. While, for PSO and GSA the elapsed time is 2013 s and 
1760 s, respectively. It is also observed that the GA approach is converged at the 28 number of iterations, 
while PSO and GSA are converged at the iteration of 38 and 75 as illustrated in Figure 4, respectively. 
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Table 2. Comparison of battery parameters based on manual, PSO, GSA and GA approaches 


Parameters Parameters extracted approaches 
(unit) Manual PSO GSA GA 
Vnom (V) 3.3000 3.2201 2.8791 2.9114 
Qnom (Ah) 1.8100 2.1164 2.1657 2.1738 
Eo (V) 3.6971 3.5246 3.6125 3.7137 
Rin (Q) 0.0165 0.0701 0.1064 0.1422 
K 0.0265 0.0038 0.0050 0.0044 
A (Ah) 0.4194 0.6806 0.6214 0.5917 
B (Ah) 4.6152 1.6488 1.4669 1.7049 
Qratea (Ah) 2.2500 
Qmax (Ah) 2.2500 
Vau (V) 4.2000 4.1351 4.1275 4.1632 
ldis,nom (A) 1.0000 
Vex (V) 3.6400 3.4287 3.3421 3.5421 
Qexp (Ah) 0.6500 1.8195 2.0451 1.7597 
Elapsed time (s) - 2013.6 1760.1 1680.1 
OF(x) - 0.1754 0.2532 0.1689 
12 
a —— PSO 
10 44 =- GSA 
5  } —S—St~S«<“*S«w ie GA 
8 41 
s i 
x | 
w 04) 
= i 
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Figure 4. The performance of the heuristic optimization methods in terms of fitness function 


In terms of accuracy, the GA approach is performing the highest accuracy with the final value of the 
fitness function, OF(x) is 0.1689, while for PSO and GSA are 0.1754 and 0.2532, respectively. Therefore, for 
further simulation study, the battery parameters obtained by the GA approach will be used. Besides that, 
Figure 5 illustrates the comparison between the discharge curve from the battery manufacturer and the 
simulated discharge curves generated by PSO, GSA and GA approach, respectively. From the figure, it is 
proven that by using the optimization method, the battery model accuracy significantly increased. 
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Figure 5. A comparison between the discharge curve from the battery manufacturer and simulated discharge 
curves generated by PSO, GSA, and GA optimization approaches 


Int J Pow Elec & Dri Syst Vol. 11, No. 1, Mar 2020 : 333 — 341 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 g 339 


4.2. Validation of BES model 

The developed Li-ion battery model utilizing optimal battery parameters obtained from the GA 
algorithm have been reconstructed to validate its discharge curves performances at different discharge 
currents. The discharge current range considered is based on nominal rating at 1C, higher (2C) and lower 
(0.2C). The results with respect to the manufacturer’s data are as shown in Figure 6. Based on the results, it is 
observed that the discharge characteristics curves of developed Li-ion model match very well with the typical 
discharge characteristic curves of Panasonic Li-ion CGR18650AF with the overall accuracy of up to 98%. It 
clearly shows that the accuracy of the developed model depends on the precision of the extracted parameters 
data from the typical discharge curves. 
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Figure 6. Validation of the discharge characteristics between the developed Li-ion battery model in 
matlab/simulink and real battery model (3.3 V, 2.25 Ah) 


5. CONCLUSION 

An easy to use and accurate battery model parameter identification strategy is introduced using a 
heuristic optimization approach. The Tremblay’s Li-ion battery model is developed in Matlab/Simulink 
based on the formulated equations, in which the parameters are obtained from different optimization 
approaches of PSO, GSA and GA algorithms. From the optimization results of the model developed based on 
Panasonic Li-ion CGR 18650AF battery test data, the GA algorithm showed the best performance with the 
fitness function and elapsed time of 0.1689 and 1680 s, respectively. The accurate battery model obtained by 
the GA algorithm was further validated by comparing the discharge curves of the developed battery model 
with the discharge curves from the manufacturer. The results showed that the GA-based optimal battery 
model gives an accuracy of up to 98%. Therefore, the obtained battery parameters by using the GA algorithm 
can be used in further simulation study related to the controller design of the lithium batteries. 
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